| CALOM 
| GV GO 
ers iri ge. 
ee. 


LOW FLOW CHARACTERISTICS 
IN ONTARIO 


MAIN REPORT 


OCTOBER 1990 


7 Environment | 
Environnement 


Ontario Jim Bradley, Minister/ministre. 


DISCLAIMER LD & 


This report was prepared for the Ontario Ministry of the 
Environment as part of a ministry-funded project. The views and 
ideas expressed in this report are those of the author and do not 
necessarily reflect the views and policies of the Ministry of the 
Environment, nor does mention of trade names or commercial 


products constitute endorsement or recommendation for use. 


Aes 
a ae 
Rae ge NRC Ri 
Perky | Lp oc me 
ee, 


' 
ra 

' 

CELA ot id 

7) a r. : ¥yT 
sh 
i 

Pare es Leip 


ry 


y, j A of t — 
at rvs ah Pa 
; A y a. 
\ * ; 1 
iy ai A x 
re . 
Meg ‘yas 
5 eh 1 
wh al j 
VW n j 
oi hey Sore 
2) i i 
: te o ‘ { Ki 
y é * ‘ 
% 
Y ‘ 
, 1 
j 
? “i> T+ ct ede “ ¥ 
ol Ae. a4 ef a as \ a, x 
pt 
fm, 4% . 
ae | Pi paw i ie ta u J be 
oad - ‘ave as a i 
5 b bi) Sad 
a 4 P “ i ~ rt 
ies oe aws int 
, a 
va 4 - ce ~ 
* i af a i eee Nene 
. 
um | 


L 


ISBN 0-7729-6827-6 


LOW FLOW CHARACTERISTICS 
IN ONTARIO 


MAIN REPORT 


Report prepared for: 
Water Resources Branch 
Environmental Services 
Ontario Ministry of the Environment 


Report prepared by: 
Cumming Cockburn Limited 
145 Sparks Ave. 
Willowdale, Ontario 
M2H 2S5 


OCTOBER 1990 


ue 
PRINTED ON 
RECYCLED PAPER 
IMPRIME SUR 
DU PAPIER RECYCLE 
Copyright: Queen’s Printer for Ontario, 1990 


This publication may be reproduced for non-commercial purposes 
with appropriate attribution 


PIBS 971E01 
log 88-2309-026 


aint ‘ 
7 wy : aT 
ve 0 ve 4 A A 
7 rh H e Te 0 Baa ; vn iD 
; a ' Ma ne ou 
a ee Te aaigentar: od Woy, Wed! 


o : ana, \ ids teen. ae. 


eye eee Pies utaR 


; dcr 

| ie | 
. . aires, if noe ear TehaSe" Sadi 
tin be a ee i eo 


x i) ae ao wn) . . | | 
ae \ ee ‘ 
2 (oe Laie hersqera Ssh 
Ler a » BOs bay’ KGET B as 
oe A wy TGQ cht 
ib si ans IHD es ‘abd CI 


"a, } v tt | eS, oad ot : i 


ACKNOWLEDGEMENTS 


The information and conclusions presented in this report were derived 
with assistance from several individuals and organizations. 


The following members of the Water Resources Branch and Environment 
Canada provided significant input and direction throughout the study: 


- Dr. L. Logan River Systems Section, Water Resources Branch, 
Ministry of Environment 

Mr. B. Whitehead River Systems Section, Water Resources Branch, 
Ministry of Environment 


- Mr. P. Pilon Hydrology Division, Water Resources Branch, 
Inland Waters Directorate 
- Mr. D. Anderson Water Survey of Canada Division, Inland Waters 


Directorate 


The background information presented in this report was obtained mainly 
from the Water Survey of Canada and the Ministry of Environment, Water 
Resources Branch. 


We would also like to express our appreciation for the time and effort of 
all others who contributed to this project by way of information, 
discussions and otherwise. 


Digitized by the Internet Archive 
in 2024 with funding from 
University of Toronto 


httos://archive.org/details/39260309040108 


LOW FLOW CHARACTERISTICS 
IN ONTARIO 


TABLE OF CONTENTS 


VOLUME 1 : 
LIST OF TABLES 


ELST. OF 


FIGURES 


MAIN REPORT : 


0 


2.0 


3.0 


4.0 


5.0 


INTRODUCTION 
1.1 General 
1.2 Study Objectives and Scope 
METHODOLOGY 
2.1 General 
2.2 Data Base 
2.3 Data Analysis and Screening 
2.4 Extreme Value Analysis 
2.5 Flow Duration Analysis 
2.6 Seasonal Analysis 
2.7 Maps 
STUDY RESULTS 
3.1 General 
3.2 Data Base Screening 
3.2.1 Data Base | 
3.2.2 Screening Results 
3.3 Extreme Value Analysis 
3.3.1 General 
3.3.2 Annual 
3.3.3 Seasonal 
3.4 Flow Duration 
3.4.1 Annual 
3.4.2 Seasonal 
3.5 Maps 
CONCLUSIONS AND RECOMMENDATIONS 
4.1 Conclusions 
4.2 Recommendations 
REFERENCES 


Page No. 


3-6 
ce! 
357 
ae 
3-8 
3-8 
4-1 
4-1 
4-2 
5-1 


APPENDIX A 


APPENDIX B 


APPENDIX C 


TABLE OF CONTENTS (cont'd) 
APPENDICES 


General Appendices 
A.l Extreme Value Analysis 
A.2 Flow Duration Analysis 
A.3 Statistical Tests - Independence, Trend, 
Randomness 
A.3.1 Test Description 
A.3.2 Example of Tests 
A.3.3 Graphical Representation 
A.3.4 Summary of Test Results 
Central Region 
B.1 General 
B.2 Station List and Record Characteristics 
B.2 Data Analysis and Screening 
B.4 Extreme Value Analysis Summaries 
B.4.1 Annual Tables 
B.4.2 Graphs of Extreme Value Analysis 
B.4.3 Summary Table of Minimum Annual and Monthly 
Consecutive 7-Day Low Flows with a 20-Year 
Recurrence Interval 
B.5 Flow Duration Analysis Summaries 
B.5.1 Annual and Monthly Flow Duration Summary Tables 
B.5.2 Annual Flow Duration Graphs 
6 Revised Analysis 
7 Maps | 
utheastern Region 
-1 General 
2 Station List and Record Characteristics 
2 Data Analysis and Screening 
4 Extreme Value Analysis Summaries 
C.4.1 Annual Tables 
C.4.2 Graphs of Extreme Value Analysis 
D.4.3 Summary Table of Minimum Annual and Monthly 
Consecutive 7-Day Low Flows with a 20-Year 
Recurrence Interval 
C.5 Flow Duration Analysis Summaries 
C.5.1 Annual and Monthly Flow Duration Summary Tables 
C.5.2 Annual Flow Duration Graphs 
Revised Analysis 
Maps 


a”QQm 
e 
N DD 


APPENDIX D 


APPENDIX E 


APPENDIX F 


TABLE OF CONTENTS (cont'd) 
APPENDICES | 


Southwestern and West Central Region 


D.5 


D.6 
D.7 


General 

Station List and Record Characteristics 

Data Analysis and Screening 

Extreme Value Analysis Summaries 

D.4.1 Annual Tables 

D.4.2 Graphs of Extreme Value Analysis 

D.4.3 Summary Table of Minimum Annual and Monthly 
Consecutive 7-Day Low Flows with a 20-Year 
Recurrence Interval 

Flow Duration Analysis Summaries 

D.5.1 Annual and Monthly Flow Duration Summary Tables 

D.5.2 Annual Flow Duration Graphs 

Revised Analysis 

Maps 


Northeastern Region 


E.5 


E.6 
E.:7. 


General 

Station List and Record Characteristics 

Data Analysis and Screening 

Extreme Value Analysis Summaries 

E.4.1 Annual Tables 

E.4.2 Graphs of Extreme Value Analysis 

E.4.3 Summary Table of Minimum Annual and Monthly 
Consecutive 7-Day Low Flows with a 20-Year 
Recurrence Interval 

Flow Duration Analysis Summaries 

E.5.1 Annual and Monthly Flow Duration Summary Tables 

E.5.2 Annual Flow Duration Graphs 

Revised Analysis 

Maps 


Northwestern Region 


General 

Station List and Record Characteristics 

Data Analysis and Screening 

Extreme Value Analysis Summaries 

F.4.1 Annual Tables 

2 Graphs of Extreme Value Analysis 

3 Summary Table of Minimum Annual and Monthly 
Consecutive 7-Day Low Flows with a 20-Year 
Recurrence Interval 


F.4. 
F.4. 


F.5 


on ee | 
e 
N ®D 


TABLE OF CONTENTS (cont'd) 
APPENDICES 


Flow Duration Analysis Summaries 

F.5.1 Annual and Monthly Flow Duration Summary Tables 
F.5.2 Annual Flow Duration Graphs 

Revised Analysis 

Maps 


TABLE 1 
2 


FIGURE 1 
fa 
3-5 
6 


LIST OF TABLES 


Stations with Unacceptable Low Flow Data Set 
Summary of Statisics for Data Used in 
Extreme Value Analysis 

Example of Extreme Value Tabular Format 

Unit Area Average Low Flows 

Example Summary of Flow Duration Analysis 


LIST OF FIGURES 


Example of Low Flow Frequency Curves 
Example of Annual Flow Duration Curve 
Example of Seasonal Flow Duration Curves 
Station Information Legend 


After 
Page 


2-2 


3-4 
3-5 
3-6 
7 


ay - ning 


an tbat sisters “ we “4 t 

at ; = a) wo: cgitevh nll h 

tee, Ween persed woh to yretabe atom = 
: a 7 ’ : . - : 


eamuola 20 Te13 


a 


at srw. yznsips74 wor Od? 26 L 3auols 
f-é See nOtsATUD WO'F Tail 0 = 
zt suind oie wetd tanod am 3 at 
GE, tig ps) eres 3 . 
o< , . 


120 


Let 


1-1 


INTRODUCTION 
General 


The knowledge of the hydrologic conditions which exist during low 
flow periods can be of primary importance in undertaking a variety 
of water resources investigations. When analysing water quality 
conditions, the low flow characteristics within the watercourse are 
a major concern to both the user and the Ontario Ministry of the 
Environment. Some uses of low flow information include the 
following: 


i) The analysis of municipal and industrial effluent discharg- 
es to streams 
ii) Instream pollutant analyses (point and non-point sources 
including mixing zones) 
iii) Reservoir design (low flow augmentation) 
iv) Environmental approvals 
v) Feasibility of small hydro developments 
vi) Water supply and evaluation for water taking permits 
vii) Base flow/groundwater recharge and/or contamination analy- 
sis . 
viii) Stream fisheries assessments 
ix) Analyse effects of changes in watershed on low flows (e.g. 
deforestation, urbanization) 
x) Agricultural 
xi) Other 


The identification of low flow characteristics within a watercourse 
is most easily accomplished using continuous hydrometric data 
recorded for the stream. Available data were previously analysed 
and published as a series of “Low Flow Characteristics" maps which 
were published several years ago (more than 10 years in most cases) 
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by the Ministry of the Environment (Ministry of the Environment, 
1973-1978). Additional data (both temporally and spatially) are 
now available and it was, therefore, considered appropriate to 
update the available data base describing low flow characteristics 
across the Province of Ontario. 


YSZ Study Objectives and Scope 


The overall objective of the proposed investigation was to carry 
out an analysis describing the low flow characteristics at suitable 
Water Survey of Canada streamflow locations in Ontario. The foll- 
owing points summarize the scope of the investigations: 


1. Update statistical low flow analyses for each of the follow- 
ing administrative regions: 
- Central 
- Southeastern 
Southwestern and West Central 
- Northeastern 
- Northwestern 
2. Undertake data base screening analyses to identify 
constraints on the data base 


3. Produce extreme value analyses for suitable stream gauging 
locations (1, 3, 7, 15, 30-day durations) with greater than 
10 years of data 

4. Produce annual flow duration analyses, using daily data for 
suitable stream gauging locations 

5. Produce monthly low flow analyses 


2.0 


2.1 


2.2 


METHODOLOGY 
General 


The available streamflow data base was obtained for all relevant 
Water Survey of Canada monitoring locations (see discussion in Sec- 
tion 2.2). Data analyses and screening were then undertaken in 
order to assess the usefulness of the data at each location prior 
to undertaking statistical analyses (see Section 2.3). Extreme 
values for selected low flow durations were then calculated for 
various recurrence intervals for both annual and monthly series 
using the techniques discussed in Section 2.4. Flow duration anal- 
yses were also undertaken utilizing the daily data base on both an 
annual and a monthly basis (see Sections 2.5 and 2.6). In most 
cases the extreme value analysis results should be used for low 
flow characteristic analysis. However, under certain conditions 
flow duration results of heavily regulated stations may give more 
conservative results. Computer drawn maps depicting low flow 
characteristics were then produced in order to summarize selected 
low flow statistics ona regional basis (see Section 2.7). 


Data Base 


The entire Water Survey of Canada daily streamflow data base for 
the Province of Ontario for the period of record to the end of 1986 
was obtained in computer tape format. 


Stations with a minimum of ten years of record were considered to 
have an appropriate record length for the purpose of this investi- 
gation. Stations inactive prior to 1981 were not considered in the 
analysis. 


Several data management programs developed by Cumming Cockburn 
Limited were then used for extraction and analysis of data from the 
computer archives. First a program was written and executed to 


Zao 


2-2 


check and screen the available data base at each station for miss- 
ing data. A second program was then run to utilize this analysis 
to compute running averages for various durations and to pick the 
corresponding minimum annual low flow. 


A third program sorted and analysed the available data to determine 
low flow values on a monthly basis. Average flows were determined 
and extracted for the annual consecutive low 1, 3, 7, 15 and 30 day 
durations and are available as part of the background files. 
Another program was used to compute flow duration curves as dis- 
cussed in Section 2.5. 


Data Analysis and Screening 


The data were manually screened to remove station records which 
contained a significant number of zero low flow occurrences for all 
durations. These stations are listed in. Table 1 and the screening 
results are discussed in Section 3.2. | 


The initial version of the Low Flow Frequency Analysis program was 
obtained from Environment Canada for use in the investigation 
(Pilon and Jackson, 1988). This program contains a set of data 
analysis and checking modules suitable for data screening and anal- 
ysis (Pilon et al, 1985). At present there is some question con- 
cerning the reliability of the methodology for checking Independ- 
ence for low flow data series (personal communication, P. Pilon, 
Water Resources Branch, Environment Canada). 


The data screening process was undertaken to test the assumption 
that the data are reliable measurements and are mutually exclusive 
and represent independent random events which are free from trend. 
The non-parametric tests applied include the Spearman Rank Order 
Correlation Coefficient for Independence; the Spearman Rank Order 
Correlation Coefficient for Trend; and a general randomness test. 
The statistical tests were applied at the 1% and 5% levels of sig- 
nificance. (See Appendix A.3 for more details on the statistical 
tests.) 
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It must also be expected that a certain number of stations which 
are tested will randomly "fail" any given statistical test. There- 
fore, the binomial distribution was used to calculate the number of 
expected random failures over the entire data base for each region 
in order to assess the suitability of the low flow data on a 
regional basis. A tabulation of theoretical and actual "failure" 
results is summarized and discussed in Appendix A.3.4. 


The effect of regulation was examined by further subdivision of the 
data base according to regulated and non-regulated stations. These 
results are tabulated in Appendices B to F and discussed in Section 
SniZs 


Extreme Value Analysis 


An extreme value analysis was undertaken for each of the 1, 3, 7, 
15 and 30-day durations for each of the stations. 


For example, the minimum annual consecutive 7 day average low flow 
was determined for each year for the station record being analysed 
(see Section 2.2). The corresponding set of consecutive 7 day 
average low flows, therefore, represents an extreme value series to 
which a theoretical extreme value distribution can be fit for the 
purpose of determining the low flows corresponding to various 
recurrence intervals. For example, the 7Q2Q low flow is often of 
interest for water quality investigation (i.e. the consecutive 7 
day average low flow with an average recurrence interval of once in 
20 years). 


Previous investigations have found that the Gumbel III distribution 
has resulted in the best fit for extreme value analysis of low 
flows on most Canadian rivers. However, for samples with large 
negative skewness, the 3-Parameter Log Normal Distribution has pro- 
ven to give adequate results and was adopted for these stations 
(Condie, 1983). A technical discussion of the fitting procedure is 
given in the users manual and briefly in Appendix A. 


2-4 


Parameter estimation for fitting of distributions proceeded in 
order of maximum likelihood, smallest observed drought and 
moments. The procedure utilized for analysis of each station 
record was identified on summary tables together with the results 
of the analyses (see Section 3.3. and Appendices B-F). 


208 Flow Duration Analysis 


The daily discharge data was extracted for each station and sorted 
by computer in order to derive empirical flow duration curves. The 
flow duration curves were summarized both numerically and graphic- 
ally for each station (See Section 3.4 and Appendices B-F). The 
tabular summary presents the results of the flow duration computer 
Output at 1% intervals across the flow duration curve. 


2.0 seasonal Analysis 


In addition to annual computations referred to above, the flow 
duration curves and extreme value analyses were undertaken ona 
monthly basis to provide an information base for seasonal water 
resources investigations. 


In order to reduce the amount of data tabulation, it was decided 
that only the results of the 7Q29 values would be summarized from 
the extreme value output. The results from the other 7-day 
duration recurrence intervals are available as part of the back- 
ground files. 


Flow duration curves were also tabulated on a monthly basis for 
comparison to the annual curves. A direct comparison with the 
plotted annual flow duration curve is, therefore, possible by using 
the numerical tabular summary. Typical results are discussed in 
Section 3.4. 


Del. 
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Maps 


Computer plotting was utilized as a cost-effective way of graphic- 
ally summarizing selected flow characteristics. Base maps were 
digitized for each region. Data overlays were then prepared in 
order to produce two maps for each region, the first map summariz- 
ing selected low flow characteristics at each station, and the 
second map providing a corresponding discharge rate per unit area. 


The resulting maps are discussed briefly in Section 3.5 and specif- 
ic maps for each region accompany Appendices B to F. 
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STUDY RESULTS 
General 


This section summarizes general study findings and refers to spec- 
ific results given in separate appendices for each of the five 
regions. (Appendices B to F). 


Data Base Screening 


Data Bases 


For all regions, a total of 389 stations were found to have 10 
years of record with the station being active within the last 5 
years. 


Other characteristics of the data base were also tabulated in Sec- 
tion 2 of Appendices B to F including, a brief station description, 
the period of record at each location, whether stream flows in the 
watershed are considered to be regulated or natural, and the drain- 
age area of the watershed. 


All of this information was extracted from the Water Survey of 
Canada HYDEX computerized data file maintained by Environment 
Canada. Statistical analyses were undertaken over the period of 
available record for both natural and regulated stations. However, 
subsequent to the data testing and extreme value analyses it was 
found that eleven data series classed as "regulated" included a 
portion of natural flow data. These stations were re-analysed with 
results addended in each Appendix as appropriate. Additional 
research indicated that three stations not recorded as mixed (nat- 
ural and regulated flows recorded for the station) in the HYDEX 
file did consist of mixed periods. An analysis of these stations; 
02GB001, Grand River below Brantford, 02GA003 Grand River at Galt, 
and, 02GA015 Speed River below Guelph, were re-analysed using only 
the regulated period of data. The period used was from 1945 (the 


3-2 


year after the Shand Dam was built) to the present for the Grand 
River Stations and 1975 to present for the Speed River data series 
(subsequent to the Guelph Dam construction). The results for the 
mixed 7Q2, 7Q19 and 7Q2q compared to the regulated period only are 
summarized as follows: 


Mixed Data Regulated aoe 
i 7 7 7. 7 
Station AP 1340 ‘390 aoe a a7 


02GB001 13.308 7.037 5.495 14.980 Sari TAS: 
02GA4003 5.822 2.176 1.679 9.186 5.120 4.156 
02GA015 1.080 05599 0.488 T2399 0.885 0.775 


The results indicate an increase in the low flow values when using 
the regulated period only, compared to the analysis of the mixed 
record period. This may indicate low flow augmentation, or trend 
in the data series. 


Users referring to analysis results for other regulated stations 
should investigate this aspect in more detail and analyses for all 
regulated stations used with care. 


3.2.2 Screening Results 


i) Zero Flows 


Prior to undertaking the extreme value analyses a total of 23 data 
sets containing a significant number (*) of zero low flows for all 
durations were identified and removed from the data base (see Table 
1). For example, for those stations where 30% or more of the 
available extreme flow data set was comprised of zero flow, it was 
arbitrarily assumed that use of this data base would result ina 
biased extreme value analysis. In addition, for short data sets 
(i.e. 10 years), reduction by 3 years results in insufficient data 
(7 years) with which to accurately fit the extreme value 
distribution. 
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A few stations which originally appeared to have 10 years of 
records were found to have a significant number of missing years in 
the detailed screening process. Those stations with minimal data 
were removed from the analysis (refer to Table 1). 


ii) Arithmetic Overflow 


Initial applicaton of the LFA program encountered an arithmetic 
overflow for a number of stations, which are identified in Table 
1. The program was subsequently modified and the analysis results 
are included as addendums to Appendices B to F where appropriate. 


111) High Outliers 


Initially the analysis of several stations with very large low 
flows (greater than 1000 m3/s) could not be completed. These sta- 
tions are identified in Table 1. Additional research was under- 
taken and the LFA program was subsequently modified to permit fitt- 
ing of the extreme value distribution. These results are addended 
where appropriate to Appendices B - F (see Table 1). 


iv) Statistical Tests 


Statistical data analysis tests were undertaken as outlined in 
Section 2.3 and described in more detail in Appendix A.3. The 
available test statistics were recently made available as part of 
the LFA (Pilon, 1988) low flow analysis program. Previous 
application of these testing procedures to low flow series has been 
limited to date. 


In general, it was found that a significant number of stations 
failed the non-parametric tests. Therefore, taken over the entire 
data base, application of these tests has indicated that the 
available data base of extreme low flows may exhibit some trend 
with some possibility of non-random characteristics. 


The data were further analysed by subdivision of the available data 
set according to length of record (i.e. >20 years and <20 years) 
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and according to regulation code. However, it was found that 
neither the length of record, nor the possible effects of regula- 
tion could account for the conclusions of the test results. One 
explanation could be that the available record lengths are too 
short to permit reasonable application and interpretation of these 
non-parametric test results. A stronger possibility is that the 
available low flow data sets do exhibit trend and non random char- 
acteristics, which could possibly be attributed to slow cyclic 
change in groundwater levels or to climatic trends. 


Additional testing was beyond the scope of the current investiga- 
tions. However, further studies are recommended since these 
results may call into question the basic assumptions underlying 
application of the extreme value analysis technique for analysis of 
low flow characteristics. 


v) Data Base 


Extreme value analyses were undertaken on an annual basis for 344 
stations and on a monthly basis for 330 stations. Further analysis 
was completed on the remaining stations with the revised program or 
by manual fitting on the annual data series. 


Fourteen stations were assigned to more than one region since they 
may be representative of either region. The fourteen stations 
which appear in more than one region are: 


Station South Southwestern North North 
die Dis Central Eastern West Central Eastern Western 


02HBO11 V 
02HB012 
02HB013 
02HB016 
O2HK002 
02HK003 
O2HK004 
O2HK005 
O2HKO06 
O4HAOO1 
04JA002 
04JC002 
04JC003 
04JG001 


SSS SSS 
SESS 
<<< 


See ee 
<< <tc e 


TABLE 2 


Summary of Statistics for Data Used in the Extreme Value Analysis 


Standard 
Deviation 


Coefficient| Minimum 


*No. of - Skew of 


Region Stations m?/s Variation 

All 82 cal 4.07 28 
78 66 Pay A =| 28 

0/5 -62 6.61 28 

78 59 UU 28 

84 MY 4 N/A 28 

Northwestern “Ths" -76 7.48 31 
we we 8.79 shi 

67 -68 9.47 31 

265 65 9.47 Sil 

abe be On 31 

Northeastern 79 -70 3.94 31 
78 64 6.25 31 

-62 58 7.68 31 

Ale: 53 8.87 31 

59 -48 10.49 31 

Southwestern otf! -67 1.24 26 
& West 74 -62 1.58 26 
Central SIH) -58 1375 26 
85 56 1.88 26 

1.03 <7 1.97 26 

Central 239 45 252 24 
42 42 Acs: 24 

244 -40 81 24 

sey 41 1,01 24 

-61 -40 1.18 24 

Southeastern 1.62 7205 12.22 28 
1.50 96 18.58 28 

1.48 93 20.97 28 

1.45 -89 22e51 28 

1.57 -90 N/A 28 


SOME STATIONS APPEAR IN MORE THAN ONE REGION. 
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The SPSS (Statistical Package Social Science, Norusis, 1986) was 
used to produce general statistics of the data base including the 
mean, standard deviation and coefficient of skew of the available 
low flow samples for different durations. These general statistics 
are summarized in Table 2. 


With reference to Table 2, it is evident that the average of the 
mean low flow increases as the duration of low flow increases. The 
standard deviation decreases as a percentage of the mean as the 
duration increases. The mean skewness of the data decreases with 
the increase in duration. 


Flow duration analyses were subsequently undertaken both on an 
annual and monthly basis. The flow duration results are summarized 
in Section 3.4, and the results tabulated on a regional basis in 
Appendices B to F. 


Extreme Value Analysis 


General 


Tables summarizing the results of the extreme value analysis for 
each region are given in Appendices B to F for various recurrence 
intervals for each duration. An example of the presentation format 
is given as Table 3. The stations are identified by the Water Sur- 
vey of Canada station number. The fitting method for the extreme 
value distribution is identified by a 3 letter code, MAX, MOM, SOD, 
PLN which stand for; method of maximum likelihood, method of mom- 
ents, method of smallest observed drought, and the three-parameter 
log normal distribution respectively. (See also Appendix A.1 for 
additional information on the fitting procedures.) 


The next three parameters appearing in the summary table are gener- 
al statistics of the data sample including the mean n-day duration 
flow (m3/s), Standard deviation, Skew (G), and the coefficient of 


TABLE 3 
EXAMPLE OF EXTREME VALUE TABULAR FORMAT 


EXTREME VALUE LOW FLOW ANALYSIS FOR 7 DAY DURATION VALUES 


208 


188 


5.8 18 20 


2.8 


RECURRENCE INTERVAL 
1.005 1.018 1.191 1.250 


STANDARD REC MIN 
C (YRS) (m3/s) 


STN§ METHOD MEAN DEVIATION G 
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398 27 
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479 1 
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Oe it ee eC 
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a ad a here 
0, 6) Neate te 
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Signe eg) ven? fe 
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er 
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No 
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7 
238 26 365 


APG Ry fe. 10 48: 


ech lb) Seriien) Le 


Wik eek ed 


C. | CME Se Fe nns 


6, bie veere 


AER ee eC 


Oy le Meera: ie. 


oF 0 ee. a) Ps 


Chet ok er 


ev 20. Tey Ho ice 


205 71 3 
923 15 
278 15 


287 68 


ey fold or © 


252 @ 
278 8 
995 @ 


Clete eae er 


1131 
43 0 
361 

696 1 


391 6 


368 2 
194111 


abde Se Ae 8. 


7 


@2KF O89 MAX 546 
@2KF@18 SOD 


O2KF OSE MAX 
@2KF 611 SOD 
@2KF@12 MAX 


oie eo Kee 


Wet [iva ne: oe 


Cag er aaa 


6, ene” oe fe 


Shee enetae ¢ 


oy Teh et Seer’ 


1.333 1.044 1.516 6.783 56 
@.218 0.106 8.811 0.485 19 
@.053 0.104 2.994 1.950 37 


@2LB005 MAX 
O2LB806 MAX 
82LB007 SOD 
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variation of record (C). The next two columns of the tables give 
the number of years of record available and the minimum observed 
average low flow for the particular duration for the available data 
set. Finally, the estimated flow for the various recurrence inter- 
ValS.arewlisted for seach «station. 


Annual 


Extreme value frequency curves were plotted for each location and 
are available in Appendices B to F. An example graph is given in 
Figure 1 which summarizes the actual data and the fitted curve for 
theel, 3557 , lo dndes0 davecurations ror the IrenteRiver at. Gien 
Ross, Water Survey of Canada Station Number O2HK004. 


A statistical summary of average unit area low flow characteristics 
was also undertaken and is summarized in Table 4. Area average low 
flow rates (litres/second/kilometre®) together with standard devia- 
tion for the 5 durations were tabulated for all the regions for the 
2 to 50 year recurrence intervals. 


It is interesting to note that in most cases, the standard devia- 
tions are approximately equal to or greater than the corresponding 
mean flows. Low flows are consistently higher than the Provincial 
average in Northwestern, Northeastern and Central Region and con- 
sistently lower in the Southwestern and West Central/Southeastern 
Regions. However, it should be recognized that this finding may be 
due to the effects of the drainage area. For example, gauged 
watersheds in Northern Ontario tend to be larger than watersheds in 
Southern Ontario. Therefore, if the relationship between watershed 
area and low flow is non-linear the corresponding statistics may be 
biased. 
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TABLE 4 


Unit Area Average Low Flows 
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Mean (Std. 


Deviation) 


Mean (Std. 
Deviation) 
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Mean (Std. 
Deviation) 
1/s/kmé 


Mean (Std. 
Deviation) 
1/s/km@ 


Mean (Std. 
Deviation) 
1/s/kmée 


1/s/kmé 


Northwestern 


Northeastern 


Southwestern 


Southeastern 


3-7 


3.3.3 Seasonal 


3.4 


3.4.1 


The 7Q2g values were also determined for each month. The resulting 
analyses are summarized in Appendices B-F. It is evident from the 
summary tables that many of the lowest of the 7Q29 generally take 
place in the summer months. 


Flow Duration 
Annual 


Flow duration tables and curves were produced to summarize the 
percentage of time the flow was greater than or equal to the given 
value. Flow duration curves for all stations are given in Appen- 
dices B to F. An example table is given in Table 5. The first 
column "Per" refers to the percentage of the period of record that 
the tabulated flow was equalled or exceeded (all flows are in 
m3/s). The annual flow duration curve values are listed in column | 
2 for the percentage summarized in column 1. Therefore, the larg- 
est daily flow recorded for this station up to 1986 is 702 m3/s, 
found in the 1st row at 0 percent. More significantly for this 
study are the low flows summarized in the later section of the 
table which have been exceeded 90, 95 and 100 percent of the time 
for the period of record. 


The 50 percent value (114 m3/s) refers to the median daily flow for 
the period of record of the flow series and can be compared to the 
mean value 145 m3/s summarized in the last row of the table. The 
annual flow duration curve is graphically depicted in Figure 2, 
which corresponds to the numerical values summarized in columns 1 
and 2. 


EXAMPLE SUMMARY OF FLOW DURATION ANALYSIS 


SUMMARY TABLE FROM FLOW OURATION ANALYSIS 


YEARS OF RECORD: 


PER 


OOoOonN @DMmMne wn —- Oo 


SBESEGESRK=>H BBSVGRBRLKBKRLB BBANBRRBRVB 


73 STATION AREA: 


AMMUAL = JANUARY = FEBRUARY 


MARCH 


TABLE 5 


O24K004 TRENT RIVER AT GLEN ROSS 


12000 km2 


APRIL 


MAY 


175.000 
173.000 
170. 000 
165 . 000 


SUMMARY TABLE FROM FLOW OURATION ANALYSIS 


YEARS OF RECORD: 


PER ANNUAL 
114.000 


111.000 
109.000 


75.300 


BRARREBR2B BSBEYBRKELKLS 
8 
8 


SS26R988 ABH 
S8dB88s8 Bas 


48. 400 


SBRIBKRERBSLK BSSRRLBKSLB dadtadaAsarnra 
2 
8 


P 
: 


JANUARY 


146.000 
143.000 
141.000 
140.000 
139.000 
138.000 
138.000 
133.000 
133.000 
131.000 


130.000 
129.000 
127.000 
127.000 
125.000 
125.000 
123.000 
122.000 
120.000 
120.000 


119.000 
119.000 
118.000 
118.000 
116.000 
115.000 
114.000 


73 STATION AREA: 


FEBRUARY 


138.000 
135.000 
134.000 
133.000 
132.000 
131.000 
130.000 
123.000 
128.000 
127.000 


154.392 


WARCH 


162.000 
158.000 
154.000 
152.000 
145.000 
142.000 
137.000 
138.000 
135.000 
132.000 


131.000 
130.000 
128.000 
128.000 
127.000 
126.000 
125.000 
123.000 
121.000 
120.000 


118.000 
117.000 
116.000 
114.000 
114.000 
112.000 
111.000 
108.000 
108.000 
108.000 


105.000 


TABLE 5 (con't) 


O2K004 TRENT RIVER AT GLEN ROSS 


12000 km? 


APRIL 


387.000 
349.000 
348.000 
340.000 
337.000 
331.000 
$26 .000 
315.000 
308 . 000 
303.000 


301.000 
734.000 
289.000 
283.000 
273.000 
285.000 
262.000 
35.000 
248.000 
245.000 


738.000 
733.000 
227 .000 
273.000 
219.000 
212.000 
204.000 
198.000 
183.000 
184.000 


182.000 


104. 000--- 478.000 


103.000 
102.000 
99.500 
97.400 
96.300 
98. 300 
$2.900 
89.500 


88. 400 
85. 100 
84.100 
81.300 


174.000 
168.000 
161.000 
157.000 
155.000 
151.000 
148.000 
142.000 


129.000 
121.000 
112.000 
108.000 


WAY 


213.764 339.733 200.326 


JUNE 


SEPTEMBER OCTOBER 
44.100 74.500 
43.300 73.800 
43.000 71.400 
42.500 69.300 
41.900 67.100 
41.00 66.000 
41.100 65.100 
40.500 63.100 
39.900 81.700 
38.500 60.300 
37.900 59.500 
37.400 58.700 
38.700 57.800 
36.500 57.200 
36.200 58.000 
3§.700 55.500 
38.400 54.600 
3.000 853.500 
34.5900 52.700 
34.00 52.100 
34.00 51.500 
33.700 50.700 
33.400 S0.100 
32.800 49.600 
32.90 48.100 
32.200 46.900 
31.700 46.200 
31.00 4§.000 
31.100 43.600 
30.900 42.200 
30.600 40.800 
279.700 39.400 
279.400 37.900 
28.600 37.400 
278.000 38.500 
27.400 33.700 
2.700 §8=633.100 
26.200 32.600 
25.700 30.900 
25.100 30.300 
24.500 29.700 
24.400 729.700 
24.100 29.400 
73.400 29.200 
73.000 28.600 
22.800 28.200 
22.500 27.700 
21.700 26.900 
21.200 25.200 
19.100 21.700 
18.200 15.000 
61.885 95.104 


34.800 


#288 £888 


Sl < 
BSBBSB BBBE 


69.800 


3.4.2 Seasonal 
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The flow duration values were also determined and tabulated for 
monthly flows. Appropriate summaries for all stations are given in 
Appendices B to F. This is also shown in Table 5 in the example 
for Station O2HK004, Trent River at Gien Ross. 


The seasonal duration curves were not graphically produced for each 
station. However an example plot for Station O2HK004 demonstrates 
how the tabulated values could be used to plot monthly curves and 
provide a comparison to the annual curve (see Figures 3, 4 and 5). 


Maps 


Selected low flow characteristics were summarized on maps for each 
region. The stations are located at the point of discharge meas- 


- urement and the selected data is summarized in an information box, 


the format of which is depicted in Figure 6. 


The stations are identified by the 7 digit Water Survey Station 
Number and this is followed by a regulation code. A code "R" indi- 
cates that data collected at the station is affected by regulation; 
the code "N" means the station data are natural or non-regulated. 
The symbols 7Q2, etc. refer to the average minimum consecutive 
7-day flow (m3/s) with a recurrence interval of 2, 5, 10 and 20 
years, followed by the minimum one day flow and the period of 
record for the station. The values shown on the right are the 
flows (m3/s) equalled or exceeded for the available period of 
record 5, 50, 75, 95, and 99 percent of the time. The symbol AREA 
refers to the station drainage area in km. 


Station names are listed along with the station numbers for identi- 
fication purposes, on the map for each region. 
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A second map summarizes the above noted low flow characteristics 
expressed as 1/s/km@. It may be possible to refer to watersheds 
with similar unit runoff rates as a means of providing preliminary 
low flow estimates for ungauged watersheds. However, the limita- 
tions of area proration should be recognized. For example, a pre- 
liminary Low Flow Regionalization Study for Southwestern and West 
Central Ontario Region (Cumming Cockburn Limited, 1988) has found 
that other significant watershed parameters enter into the deter- 
mination of low flows. Additional investigations are required in 
order to refine estimating techniques for ungauged watersheds. 


4-1 


4.0 CONCLUSIONS AND RECOMMENDATIONS 


4.1 Conclusions 


1. 


The analysis of non-parametric test results indicate that the 
available data base of extreme low flows may exhibit some 
trend and dependence with some possibility of non-random char- 
acteristics. Therefore, some degree of caution should be used 
when applying the results of this study. 


The Gumbel Extreme Value Distribution was found to adequately 
fit the majority of available low flow series for various low 
flow durations. However, for a number of samples with large 
negative skewness, the 3 PLN Distribution was adopted. 


Flow duration analyses were successfully undertaken both on an 
annual and monthly basis. 


Extreme.value analysis were undertaken on an annual basis for 
344 stations and on a monthly basis for 330 stations using the 
LFA program. Graphical analyses were undertaken for the 
remainder of the stations. 


It was found that area average low flows are consistently 
higher than the Provincial average in the Northwestern, 
Northeastern and Central Regions and consistently lower in the 
Southwestern/West Central and Southeastern Regions. 


A map summarizing the following low flow characteristics was 
produced for each region: 


- 7 day average extreme values for the 2, 5, 10 and 20 year 
recurrence intervals 

- flow duration values which were equalled or exceeded over 
the available period 5, 50, 75, 95 and 99 percent of the 
time 


Ae? 


A second map was also produced for each region which summariz- 
ed the above low flow characteristics expressed as 1/s/km2. 


Users referring to the analysis results for regulated stations 
should investigate the effects of regulation on low flows in 
more detail. 


4.2 Recommendations 


1. 


Further investigation should be undertaken to confirm the app- 
licability of available non-parametric tests for low flow 
series. 


The possibility of effects of cyclic changes in groundwater 
regime or climatic changes on low flows should be examined in 
future investigations. 


Additional investigations are required in order to refine low | 
flow estimating techniques for ungauged watersheds for each 
region. 


Maps summarizing monthly 7Q29 low flow characteristics should 
be produced for each region. This would be useful when under- 
taking seasonal low flow analyses to facilitate seasonal anal- 
yses. 


Discharge data are collected continuously by the Water Survey 
of Canada at each station. Data analysis and management tech- 
nique are now available which would allow efficient updating 
of the present analyses on a frequent basis. In our opinion, 
the low flow analyses should be updated every three years in 
order to provide reasonably accurate information for investi- 
gations requiring low flow information. 
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APPENDIX A.1 
EXTREME VALUE ANALYSIS 


Introduction 


Statistical distributions allow estimates of probability of. exceedance of 
events to be made by analytical techniques. Various methods are avail- 
able, and a good discussion is given by McMahon and Mein (ref 9). 


At a gauged location, average low flows can be determined for selected 
durations of time for each year of record. For example, the corresponding 
minimum average consecutive 15-day low flow can be determined for each 
annual record. The series of 15-day average annual low flows can then be 
ranked lowest to highest and an extreme value analysis undertaken using a 
theoretical distribution. 


The Gumbel Extreme Value Distribution is commonly used for fitting low | 
flow frequency curves to the available data. 


The Gumbel Type III Distribution 


This is a variation of Fisher and Tippett's third asymnptotic distribution 
of extreme values and is sometimes referred to as the Weibull distribu- 
tion. In view of his many contributions it is often referrred to in 
hydrology as the Gumbell II. The probability density function is: 


Neate a-] ¥ ae | a (1) 
GER eee Gb fet eet 


where e is the lower bound parameter 
u is the characteristic drought 
and ais the shape parameter 


Various methods are available for determining the distribution parameters 
for a particular data set. See references 2 and 6 for example. 


The Probability Function 


The density function is integrable and gives the distribution function 


a (2) 
F (x) = 1 - exp . | Kvene | | 


which gives the probability of non-exceedance of x. 
Since it iS more common to require x for a given probability of 
non-exceedance, a simple re-arrangement gives 


1/a (3) 
X =e + (u - e) | - In bd - F(x) | } 


There are occasions when the sample series of low flows can have a large 
negative coefficient of skewness, and since the Gumbel Type III distribu- 
tion cannot have a skewness of less than -1.14, then the distribution 
cannot be fitted. | 


There are insuffucient natural samples available with such low skewness to 
make a firm choice of an alternative treatment, but from the few avail- 
able, it was found that the negatively skewed three-parameter lognormal 
distribution provided an acceptable alternative. 


If y = 1n (a-x) is normally distributed, then the probability density 
function of x is given by 


‘ 2 
x) -——_— | <5 [in (a-x) -m | (4) 
(a-x)o /21 


where m and o are respectively the mean and standard deviation of the 
series 1n (a-x) and a is an lower boundary parameter. This form of the 
distribution can only have negative skewness. The equation in the form 
In (x-a) is for distributions with positive skewness. 


Taking moments, re-arranging and replacing them by their sample estimates 
gives: 


k3 + 3k - gy = 0 (36a) (5) 
and after solving for k 


Xmas Kk (37a) (6) 


a = 
m = in(-s/k) - 1/21n(k2+1) (38a) (7) 
o2 = 1n(k2+1) (39a) (8) 


Thus, the distribution is completely defined and the T-year low flows can 
be computed from: 


QT = eee exp(m ot) (9) 


Additional details are discussed in the user's manual. 


Source: Condie, R., L. C. Cheng 
"Low Flow Frequency Analysis - Program LOFLOW", 
Water Resources Branch, Inland Waters Directorate 
Environment Canada, Ottawa, 1987 
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FLOW DURATION ANALYSIS 


APPENDIX A.2 
FLOW DURATION ANALYSIS 


A flow duration curve is a plot of the flow of a stream against the per- 
cent of time the indicated flow was equalled or exceeded during the period 
covered by the available flow data. This curve is extremely useful for 
hydro power studies and for characterizing the local streamflow regime. 


An empirical procedure was used to analyse each sample discharge record. 
Flow duration curves are derived by rearranging the available daily or 
monthly flow data in order of magnitude. The total time period represents 
100% of the time. Therefore, by definition of the procedure used here, 
the largest value is exceeded 0 percent of the time and the smallest value 
is exceeded 100% of the time. 


The flow duration curve represents the percent of time that a specific 
discharge occurs at that location. However, the curve does not indicate 
the period of time in the year when the flow is less than or equal to the 
selected value. Therefore, in some instances, it is also useful to devel- 
op flow duration curves on a monthly basis, that is, all data for the 
month of January over the entire period of record is analysed independent- 
ly to produce a flow duration curve representative of flow conditions in 
January. 


Single station flow duration analyses and corresponding monthly flow dura- 
tion curves have been determined at all relevant stream and flow locations 
across the Province. 


; bide: fo diatone fa Ths e f es a ae 
N ees a ee oye Se ie Ch 
wi 1 oe ce 
i. toe ant to Le ae td 


a i iA 
hs 1) x Yate Figen oe | "Sah We va seis 


1 
‘ai 


i Mee i _ in ; ah do vd 
(taht saed os ‘alia ae yi waengen *0 Hie: py wat), el ftoem 
Ce ee te ee ed 
Fine teat! vat 
"i he que e papsenis zi 


ay wig all, af ph Wetigniticer!? ant 
“Vi sot siee abet ‘wie GOxbrisei bi 


& 


hdl i it anh PD: WOEsen, 400 


v4 net ¥ «WV by (hets9 | 50 
Po ‘ yids Whe eb wory an 
at wee woaunnll to stron 
wen s snow oY y 

| ‘a rainvet 

1 AOizat: stent? 

iol dawd?no!? 

SN ZzqnD« 


APPENDIX A.3 
NON-PARAMETRIC TESTS FOR 
INDEPENDENCE, TREND AND RANDOMNESS 


A.3.1 


A.3.1.1 


APPENDIX A.3 
NONPARAMETRIC TESTS FOR INDEPENDENCE, 
TREND AND RANDOMNESS 


Test Description 


This appendix briefly summarizes the functions evaluated in the 
package and gives the methods used to determine their statistical 
Significance. Statistical tables are provided for ease of 
reference. 


Introduction 


Any statistical test of significance will generally be made using 
the following steps: 


a) 


b) 


c) 


d) 


e) 


F) 


State the null hyposthesis, Ho. For instance in split sample 
tests, the null hypothesis may be that there is no difference 
between the sample means. 


Choose a significance level, a. 


Choose an appropriate statistical test. In this program all 
tests are nonparametric. 


Compute the test statistic. 
The sampling distribution of the test statistic is known and 
has been tabulated, and the chosen significance level then 


defines the region of rejection. 


If the computed test statistic lies in the region of 
rejection, then the null hypothesis is rejected. 


Consider now the four tests in this program. 


A.3.1.2 The Spearman Rank Order Serial Correlation Coefficient for 
Independence 


If the series Q; with i ranging from 1 to N is put in 
chronological order, two time series are formed and their 


respective ranks computed: 


Q1, QD ecccccccccceeeQu-1 by Xi the rank of Q:3 j=1; N-1 
ANd QD, QZ eccccccccccee Qu by Yas the rank of Q;3 i=2, N 


then Spearman rank order serial correlation coefficient is 


S)= 4 (Ex%+ Zye = Zd¢) (Ex¢ sy2)~* (1) 
where Exé = (m3-m)/12- =Tx 
By = (m3-m)/12- BT 


qd, is the difference in rank between x and Y; 
m = N-1 


and the summations are over the m pairs of Xis Vas 


Ignoring for the moment the terms in T and putting them at zero, 
equation (a) becomes 


Si 1- (62d?) (m3-m) (2) 


the more familiar form of the Spearman rank correlation 
coefficient. 


The terms in T adjust for tied ranks and are computed as 

follows. If for instance three observations in the x series were 
tied for ranks 17, 18, and 19 then each observation is given the 
rank 18; if two were tied for ranks 24 and 25, then each is 
ranked 24.5. 7 


A.3al.3 


For each tied set, T is computed from 


barsaue 
T, = (d9 - g)/12 


where J is the number of observations tied at a give rank. 
21, and 21, are defined by the extension of the foregoing. 
When N is 10 or greater, then the function 


[(m-2)/(1-S2)]? (3) 


is distributed like Student's t with m-2 degrees of freedom. A 
one-tail test must be used. 


The Spearman Rank Order Correlation Coefficient Test for Trend 


If the seriés Q, with i ranging from 1 to N is put in 


chronological order, ranks assigned and denoting the series 


Q1, 2 sevcccceeeeQn bY y,» the rank of Q; 
ANG Bh I ainae aciae we NODY Xi the sequential order of Q; 


then the Spearman rank order correlation coefficient ae is 
calculated as in equation 1, except that m=N, (i = 0, and the 


summations are taken over the N pairs of Xia Vie 


For N 


10 or greater, then the function 


ct 
u 


eke [(N-2)/(1-r2) * (4) 


is distributed like Student's t with N-2 degrees of freedom. The 
null hypothesis is that there is no trend, either upward or 
downward with time, and so a two-tail test is used. 


A.3.1.4 Runs Above and Below the Median for General Randomness 


This randomness test is based on the order or sequence in which 
the individual scores or observations were obtained. Actually, 
the test is based on the number of runs which a sample exhibits. 


A run is defined as a succession of identical symbols which are 
followed and preceded by different symbols or by no symbols at 
all. 


The total number of runs in a sample of any given size gives an 
indication of whether or not the sample is random. If very few 
runs occur, a time trend or some bunching due to lack of indepen- 
dence is suggested. If a great many runs occur, systematic 
short-period cyclical fluctuations seem to be influencing the 
sample. 


For example, once the median of the sample has been determined, 
each observation can be labelled as being above and equal to or 
below and equal to the median. If "A" represents above and equal 
to the median and "B" represents below and equal to the median, 
then a sample may be ordered as 


AABBBABBBBAABA 


(A run represents a succession of identical symbols). For our 
example, each run is underscored and numbered consecutively: 


BBB AA 


AA 88B A BBBB AA OB 
1 3 6 


ne 


This sample begins with 2 observations above or equal to the 
median, followed by a run of 3 observations below or equal to the 
median, etc. 


Seven runs are observed in all: that is, the total number of 
runs above and below the median RUNAB, is 7. If n] represents 


A.3.2 


Source: 


the number of events of type A, then ny = 6. If ng denotes the 
number below the median, type B, then no = 8. Thus, the number 
of observations is equal to (nz + N92). 


In order to apply this run test, one must determine nj, n2 and 
RUNAB. 


The null hypothesis, Ho» is that the A's and B's occur in random 
Order. The alternate hypothesis, Hs is that the order of the 


A's and B's deviates from randomness. 


When either ny or ng is greater than 20, the sampling distribu- 
tion of RUNAB tends to normality with 


|RUNAB - [(2nyn2)/(ny+n2)+1] | 
Le (5) 


{2nyng(2nyng-ny-ng)/E(ny+ng)(ny+ng-1) 1? 


where z is an N(0,1) variate as described in Table A.4. This 
package uses a region of rejection defined by 

z greater than 1.96 for a = .05 

z greater than 2.575 for a = .01 


Example of Tests 


The non-parametric test referred to in Section 2.3 were performed 
On all extreme value data sets. The data for the 7-day duration 
for Station 02BD002 which "failed" the non-parametric tests is 
tabulated in Table 1 and comparatively the data for Station 
O02AB009 which passed all the non-parametric tests is given in 
Table 2. The results of the tests for independence, trend and 
randomness for both stations are given in Tables 3 and 4 respect- 
ively. 


Pilon, P.J.,.R. Condie, K. OB: Harvey 

"Consolidated Frequency Analysis Program - CFA", July 1985 
Water Resources Branch, Inland Waters Directorate, 
Environment Canada - Ottawa, 1985 


TABLE 1 


7 DAY DURATION FLOWS USED 
FOR EXTREME VALUE ANALYSIS 


SHEBANDOWAM RIVER AT SUNSHINE O2AB009 
’ DAY LOM FLOW MEAN DISCH. IN PERIOD JAN 01 DEC 31 

STARTING ? DAY ASCENDING CUMULAT. - RETURN 
MONTH YEAR MEAN FLOU ORDER RANK PROGABIL. PERIOD 
(CNS) (CAS) ee) (YEARS) 
8 1958 9.2310 2.0900 ih 2.95 48.67 
1 1959 ?. 3610 2.2740 2 3.48 LOaco 
8 1960 2.9910 2.3940 3 8.90 bees 
8 1961 7.7160 2.6590 4 12033 S211 
10 1962 7.7990 2.3640 2 1B (6) S230 

10 1963 2.9100 2.8910 6 19.18 J.2 
1 1964 2.2740 3.2360 ? 22. 60 4. 42 
8 1965 4.1800 3. 3830 3 26.03 3.84 
? 1966 2.5940 3.9100 9 29.45 3.40 
10 1967 3. 6090 3. 6090 10 32. 88 3.04 
if 1968 3. 4610 3.7160 i 36. 30 Lets 
3 1969 4.6960 4.0070 12 39.73 LERY: 
3 1970 2.84646 4.1800 8) 43.15 fi BY 
3 1971 3. 7980 4. 6560 14 46.58 cus 
9 1372 9. 7010 3.2310 ioe 90.00 2.00 
1 1973 6.1360 3. 4270 16 93. 42 Tobe 
8 1974 §.1490 . 9.4610 17 96.89 1576 
11 1975 &. 3770 ~ 9.9010 18 60.27 1.66 
12 1976 2.0900 5. 9570 1? 63.70 Hyg 
2 1977 2.6590 9. 9966 20 67.12 1.49 
12 1978 7.7360 6.1360 at 70.55 te a2 
9 1979 4.0070 6.1490 ey (anee ass 
6 1980 3. 3830 6.5970 23 77. 40 fr? 
9 1981 3.2360 6. 7060 24 80.82 1324 
2 1982 3. 4270 7.7360 a3 84.25 Liga sits 
4 1983 10. 8860 8.2340 26 87.67 1.14 
10 1984 §. 7060 9. 3610 27 SLeL0 eo, 
L 1985 8.2340 9.7990 28 94.92 1.06 
12 1986 3.9970 10. 3860 2? o0 095 d2u2 


TABLE 2 


7 DAY DURATION FLOWS USED 
FOR EXTREME VALUE ANALYSIS 


APCHEPICOTEN RIVER AT HIGH FALLS U2EDNNZ 
? DAY LOW FLOW MEAN DISCH. IN PERIOD JAN O1 DEC 31 
STARTING ? DAY ASTERDING CURULAT. RETURN 
MONTH YEAR MEAN FLOW ORDER RANK PROBABIL. PERIOO 
(CRS) (CAS) (2) CYESRS » 
10 1924 7.0000 3. 3860 i FasL) ee 38.87 
? 1925 20.9710 9. 6090 2 2.70 37.00 
4 1926 17.0140 15. 2710 3 4. 39 iS oyatit 
? 1927 23. 6430 19. 9710 4 6.98 16 44 
3 1928 29. 7000 17. 0000 2 Peed 12.387 
9 1933 19.2710 19.0140 6 9.46 10.97 
1 1934 21. 3140 20. 4430 ? 11225 bart 
? 1935 27.3710 20.9710 8 12.24 ua 
12 1936 29. 3140 21.1570 9 14.93 6.38 
4 1937 23.1900 21.2140 10 16.22 &, 1¢ 
10 1938 23.0570 21. 3140 a 17.34 2. 08 
8 1939 23.1860 22-0010 12 19.9? 9. 10 
3 1740 25. 8000 23.0570 i3 21. 28 4.70 
8 1941 24. 2430 23.1000 14 2e.90 4.35 
rd 1942 27.9970 23.5000 15 24. 66 4.05 
12 1943 25. 9430 23.8860 16 26.35 ae 
1 1944 23.8860 24. 2430 17 28.04 7.57 
3 1945 24.5430 24.5430 18 29.73 3. 36 
10 1946 26.1570 25.1860 19 gu 3506 
12 1947 21.1570 29. 3140 20 35,40 3.92 
10 1948 29. 6430 25. 6430 21 34. 80 2.87 
12 1949 23,5000 29;-80000R") 22 36. 49 Zz. .08 
2 1950 31.1710 25. 9430 23 38.18 2.62 
3 1991 90.0860 26.1970 24 39. 86 cas 
il 1952 44.0570 26. 7000 25 41.55 endl 
10 1953 31.1900 27. 2000 26 43.24 ars 
il 1994 36.2140 27. 3710 27 44.93 rae 
9 1959 21. 2140 27.9570 28 46.62 2.14 
1996 32. 2290 27. 7000 29 48. 31 2.07 
7 1997 3?. 3710 29. 7000 30 90.00 2.00 
8 1958 31. 8900 30. 4430 31 51. 6? Aes 
6 1959 42. 4000 30. 8710 32 33. 38 1.87 
1) 1960 30. 4430 21.1000 33 95.07 1. 82 
? 1961 41.0900 71.1710 34 58.76 1. 764 
8 1962 42.1860 31.8000 35 98.45 71 
11 1963 38.2290 32.2290 36 60.14 1. 66 
6 1964 30. 8710 32. 4570 37 61.82 1.62 
9 1969 61.0140 34. 7000 38 63.91 As. 52 
10 1966 36. 6290 36. 2140 33 69. 20 1203 
12 1967 48. 7430 76. 4290 40 86.89 1. 49 
3 1968 32. 4570 38. 2290 41 68.98 1.46 
9 196? 2. 3860 39. 3710 42 70.27 1. 42 
12 1370 91.0430 39, 4430 43 71.96 1.39 
il 1971 34. 9000 3%. 7570 44 73.69 1. 36 
9 1972 45.1860 41.0000 45 79. 34 3b. 33 
LZ 1973 08. 3000 41. 4710 46 77.03 1. 30 
bs) 1974 39.7570 42.1860 47 78.72 27 
11 1975 15.9710 42. 4000 48 80. 41 1.24 
10 1976 22.7710 44.0570 49 82.09 iene 
1 1977 9. 6090 45.1860 30 83.78 £519 
10 1978 41.4710 46.4710 “pi 85.47 1.42 
? 1979 93.1430 48. 7430 32 87.16 £15 
6 1980 39. 6430 30. 0860 33 $8.85 Lek3 
3 1981 27. 7000 90. 2290 34 90.54 1.10 
N) 1982 26. 7000 91.0430 99 Fay € 1.08 
9 1983 99.1430 93. 1430 36 V3.3 1.06 
10 1984 90.2290 99. 1430 37 95. 61 1205 
s) 1985 44.4710 98. 3000 38 97. 30 1.03 
8 1986 27. 2000 61.0140 a 98.99 1.01 


TABLE 3 


RESULTS OF STATISTICAL TESTS 
FOR STATION O2AB009 


--~ SPEARMAN TEST FOR IMDEPENDERCE --- 


G2AE0091 3007 SHEEANDOWAN RIVER AT SUNSHINE 
ANNUAL MINIMUM DAILY FLOW SERIES 98 TO 86 DRAINAGE AREA = 2800. 000 


SPEARNAN RAXK ORDER SERTAL CORRELATION COEFF 
CORRESPONDS TO STUDENTS T 
CRITICAL T VALUE AT e LEVEL 


O2068 meDSF.= “26 
0.347 
1.706 NOT SICHIFICANT 
2.4977 NOT SIGNIFICANT 


te ee 


Interpretation: The null hypothesis is that the serial(lag-one) correlation 
iS zero. 


At the 5/ level of significance, the correlation is not significantly 
different fron zero. That is. the data do not display significant serial 
dependence. 


--- SPEARMAN TEST FOR TREND --- 


02A£00913007 SHEBANDOWAR RIVER AT SUNSHINE 
ANNUAL MIKIMUN DAILY FLOW SERIES 58 TO 86 DRAINAGE AREA = 2800. 000 


v 


SPEARMAN RANK ORDER CORRELATION COEFF 
CORRESPORXDS TO STUDENTS T 
CRITICAL T VALUE AT - LEVEL 


=QR227 DF saeesel 
-1.209 

-2.052 NOT SIGNIFICANT 
=25001 NOT SIGHIFICART 


Interpretation: The null hypothesis is that the serial(lag-one) correlation 
is zero. 


At the 57 level of significance, the correlation is not significantly 
different fron zero. That is, the data do not display Sigoificaat trend. 


~~~ RUN TEST FOR CENERAL RANDOMNESS --- 


02A600913007 SHEBANDOWAN RIVER AT SUNSHINE 
ANNUAL MINIMUM DAILY FLOW SERIES 58 TO 86 DRAINAGE AREA = 2800. 000 


THE RUMNBER OF RUXS ABDVE AND BELOW THE MEDIANCRUNAR) = 11 
THE NUMBER OF OBSERVATIONS ABOVE THE MNEDIANCHL) = 14 
THE NUMBER OF OBSERVATIONS GELOW THE MEDIAN(H2) = 14 
Range at 57 level of significance: 10. to 20. HOT SIGNIFICART 


Interpretation: The navll hypothesis is that the data are randon. 


At the 57 level of cp pdb the null hypothesis cannot be 
rejected That is, the sanple is significantly randon. 


TABLE 4 


RESULTS OF STATISTICAL TESTS 
FOR STATION O2AB009 


--- SPEARMNAN TEST FUR INDEPENDENCE --- 


0260002130607 MICHIPICOTER RIVER AT HIGH FALLS 
ANNUAL MIXIMUN DAILY FLOW SERIES 24 TO) 86 DRAINAGE AREA = 5130. 000 


SPEARHAN RANK ORDER SERIAL CORRELATION COEFF 
CORRESPONOS TO STUDENTS T 
CRITICAL T VALUE AT cH LEVEL 


“uu ou 

Oe) 
co 
w 
< 


Interpretation: The null hypothesis is that the serial(lag-one) correlation 
is zero. 


At the 1/ level of significance, the correlation is significantly 
different from zero. That is. the data display highly significant 
Serial dependeace. 


--- SPEARMAN TEST FOR TREXD --- 
026000213007 MICHIPICOTER RIVER AT HIGH FALLS 
ANNUAL MINIMUM DAILY FLOW SERIES 24 T0 86 DRAINAGE AREA = 5130. 000 


SPEARMAR RANK ORDER CORRELATION COEFF =0.528 “OD. F 2" 37 


CORRESPONDS TO STUDENTS T : -4. 693 
CRITICAL T VALUE AT 527 LEVEL = -2.003 SIGRIFICAXT 
= eM YAP ee EET AW Y SIGRIFICART 


Interpretation: The null hypothesis is that the serial(lag-one) correlation 
is zero. 


At the 127 level of significance, the correlation is cic 
differeat from zero. That is the data display highly significant trend 


--- RUN TEST FOR GEXERAL RANDOMNESS --- 


626000213007 NICHIPICOTER RIVER AT HIGH FALLS 
ANNUAL MINIMUM DAILY FLOW SERIES 24 10 86 DRAINAGE AREA = 5130. 000 


THE. NUMBER OF RUNS ABOVE AND BELOW THE MEDIANCRUNAB) = 11 
THE NUMBER OF OBSERVATIONS ABOVE THE MEDIANCHL) = 29 
THE NUMBER OF OBSERVATIONS BELOW THE MEDIAK(N2) = 29 
(MOTE: 2 1S THE STANDARD HORMAL VARIATE. > 
For this test, 2 = 5.034 
Critical Z value at the 57 level = 1.960 SIGNIFICANT 
Critical Z value at the 17 level = 2.575 SIGHIFICAXT 


Interpretation: The null hypothesis is that the data are randon. 


At the 12 level of significance, the null hypothesis can be rejected. 
That is, the Sanple is not significantly raadon. 


A.3.3 


A.3.4 


Graphical Representation 


Graphical procedures are available within the LFA (Pilon, 
Jackson, 1987) program which visually depict some of the results 
of the non-parametric tests as shown in figures 1 to 6. 


Figures 1 and 2 show that there is an apparent trend for Station 
02BD002 when compared with the rank versus time graph of Figure 
2. These results correspond to the statistical test results. 

The increase in low flows as a function of time is very apparent 
for the data for Station O02BD002 (see Figure 3) when compared to 
the data from Station 02AB009 shown in Figure 4. Figures 5 and 6 
graphically depict the probability density function by displaying 
a histogram of the number of occurrences of between ranges of low 
flows. It can be seen that both stations are positively skewed 
for the 7-day duration data. 


Summary of Tests Results 


Test results are summarized by region for significant levels of 
1% and 5% in appendices B to F, Section 2.0. It can be seen 
from the summary tables (Tables 2 to 6) that the data for 
stations in the Central Region are random at the 1% level (refer 
to Table 2). 


All other test results indicate that a large number of the data 
sets "fail" the non-parametric test. This could indicate some 
dependence, trend and non-randomness. Hence, some degree of 
caution should be used when applying the results of this study. 


To examine effects of length of record and regulation, the 
analysis was completed using criteria of >20 years of record, <20 
years of record and regulated vs. non-regulated station data. 
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TABLE 2 
Summary of Data Screening All Stations 


Northeastern Region 


Independence 


Randomness 


Randomness 
Day 
Duration 


30 71 27 
“ome [oe | 


Central Region 


Independence Randomness 
Day 
Duration 


Randomness 


Day 
Duration 


* Total of the 5 durations for stations in this region 


Our : The duration the data set represents {je average 30 day low flow 

Sig : The number of stations which show significant dependence, trend, non randomness 

Not : The number of stations which show independence, free from trend, and randomness 

Per : The percent binomial probability that this number of stations would fail the non parametric tests 


Tables 3 and 4 show the results of data from > 20 years of record 
with Table 3 showing less non-regulated stations than the 
regulated stations on Table 4. The regulated stations show a 
greater percentage of values which can be said to "pass" the 
non-parametric test compared to those which "failed". Although 
the Central Region data is showing computed randomness at the 1% 
level of significance, the binomial probability is. still "0" for 
all other tests. 


Tables 5 and 6 summarize the analysis results for stations with 
less than 20 years of record. The findings suggest that neither 
the length of record or the effect of regulation can account for 
the overall conclusion that the statistical tests indicate some 
degree of dependence, trend and non-randomness of the low flow 
data sets. 


It is possible that application and interpretation of the 
non-parametric tests at the 10% level of significance could 
result in fewer statons "failing" the non-parametric tests. 


TABLE 3 


Summary of Data Screening 
Non Requlated Stations With A Period of Record Greater Or Equal to 20 Years 


Northeastern Region 


Independence Randomness 


Day 


Ouration Per Sig 


10 2 


0 


11 
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Northwestern Region 


Independence 
Day Randomness 
Duration 


Randomness 


Day 
Duration 


% 
1 [ 
3 L 
7 19 
15 19 17 13 13 20 22 
30 14 10 16 20 20 
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Central Region 


Independence Randomness 
Day 5% 
Duration 


17 

17 

18 

19 
2 


Day 
Ouration 


* Total of the 5 durations for stations in this region 


Dur : The duration the data set represents te average 30 day low flow 

Sig : The number of stations which show significant dependence, trend, non randomness 

Not : The number of stations which show independence, free from trend, and randomness 

Per : The percent binomial probability that this number of stations would fail the non parametric tests 


TABLE 4 


Summary of Data Screening 
Regulated Stations With A Period Of Record Greater Or Equal To 20 Years 


Northeastern Region 


Independence Randomness 


Day 
Duration 


Northwestern Region 


Independence Randomness 


Day 
Duration 


Day 
Duration 


1 ns 0 0 
3 29 19 0) 22 0 16 32 
7 29 19 0 24 0 19 29 
2 oe es 0 0 i * 
0 0 
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Day 
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Randomness 


Day 
Duration 


0 7 els 
0 8 15 0 
0 14 5 17 a 0 
15 0 e 7 7" : 
30 0 6 
a a a 


* Total of the 5 durations for stations in this region 


NOOnw 


Our : The duration the data set represents ie average 30 day low flow 
Sig : The number of stations which show significant dependence, trend, non randomness 
Not : The number of stations which show independence, free from trend, and randomness 


Per : The percent binomial probability that this number of stations would fail the non parametric tests 


TABLE 5 


Summary of Data Screening 
Non Requlated Stations With A Period Of Record Less Than 20 Years 


Northwestern Region 


Independence Randomness 


Day 
Duration 


1 


5 
ier fee 


Randomness 
Day 
Duration 


14 0 0 14 
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Randomness 
Day 
Duration 


Central Region 


Independence Randomness 
Day % 
Duration 


Day 
Duration 


8 


* Total of the 5 durations for stations in this region 


Dur : The duration the data set represents ie average 30 day low flow 

Sig : The number of stations which show significant dependence, trend, non randomness 

Not : The number of stations which show independence, free from trend, and randomness 

The percent binomial probability that this number of stations would fail the non parametric tests 


TABLE 6 


Summary of Data Screening 
Requiated Stations With A Period Of Record Less Than 20 Years 


Northwestern Region 


Independence Randomness 


Day 
Duration 


Northeastern Region 


Independence Randomness 


3 0 6 
3 0 6 
3 0 6 
3 3 6 
3 6 
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Duration 
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Central Region 


Independence Randomness 
Day 
Duration 


Day 
Duration 


1 
3 
7 
15 


30 
* TOTAL 


* Total of the 5 durations for stations in this region 


Dur : The duration the data set represents ie average 30 day low flow 
Sig : The number of stations which show significant dependence, trend, non randomness {| 
Not : The number of stations which show independence, free from trend, and randomness . | 


Per : The percent binomial probability that this number of stations would fail the non parametric tests 


